Correction: Impact of elevated air temperature and drought on pollen characteristics of major agricultural grass species

[This corrects the article DOI: 10.1371/journal.pone.0248759.].


Introduction
Grass pollen is the major cause of aeroallergen-induced respiratory diseases [1,2]. Besides hay fever, grass pollen can also lead to severe asthmatic reactions in the lower human airways when the pollen structure is decomposed and granules are released [3].
In the last decades the length of the pollen season worldwide has been extended and just comparably, the time period in which pollen allergies occur has been prolonged [4]. These clear changes in the pollen season can be explained by an earlier flowering of plant species due to global climate warming [5,6] whereas the end of flowering seems to be largely unchanged [7]. Besides longer pollen flight seasons, the pollen production and allergenicity of pollen are influenced by higher atmospheric CO 2 concentrations [8] and an increasing number of droughts and elevated temperature [9,10]. The affected areas by such extreme conditions will increase in the next decades as well [11]. A change in the growing conditions and thus competition as well as newly invading species will also lead to a shift in the composition of species [12]. In agricultural grassland systems, extensively planted cultivars such as the early flowering variety Ivana from the species perennial ryegrass (Lolium perenne L.) are known to have a low tolerance level against drought due to their e.g. alpine origin with high precipitation quantities and thus their cultivation will be limited in the future. Thus, a wise selection of species and cultivars which are adapted to the altered climate conditions is necessary in order to maintain agricultural yields. It was shown that C 3 grasses in general have higher requirements regarding water availability than C 4 grasses, thus it is suggested that C 4 grasses would be more abundant when water limitation increases [13]. The change in species/cultivar composition and management such as cutting dates on agriculture land and other grassland types might also affect the pollen production and their allergenic potential [14]. Unfavorable conditions during the growing season such as extreme dry periods and heat typically reduce grass growth and might inhibit pollen production, while allergenicity potentially increases due to plant stress [15,16]. A study from Switzerland [17] showed that under the hot and dry conditions in spring and summer 2003 peaks in pollen concentration were already reached in May or beginning of June while the duration of the grass pollen season tended to be shorter than in other years. It was also found that grasses almost stopped growth and pollen production at an early stage by end of June. Accordingly patients allergic to grass pollen had severe symptoms of hay fever in May and June whereas symptoms were reduced towards the end of June. Temperature and precipitation are the main drivers of plant growth and pollen development [18,19]. It was already shown that herbaceous taxa such as grasses are highly climate sensitive, especially for water availability, compared to other taxa [20]. A long term study [21] on grass pollination at the western Mediterranean coast revealed that elevated minimum temperatures and a rise of precipitation in spring led to higher average pollen concentrations and an earlier ending of pollen season. Other studies reported that an increase in temperature and precipitation intensifies the pollen production of early flowering species, while there is only a small effect on late flowering species [20,22], likely related to differentially impacting cutting dates [14]. The timing of water availability during the growing season also plays an important role for the plant development and affects the number of inflorescences, as it could be shown for tallgrass [23]. Nevertheless the response to the water availability during the growing season was still species-specific [23].
Under extreme growing conditions, e.g. longer drought or warm periods, plants suffer from water stress. It has already been shown that there is a clear link between plant growth and-stress of mesic temperate grasslands [13]. Whether the grass allergens are influenced by plant stress remains unclear, since their basic function inside pollen is unknown for the majority of those allergens. According to other studies the pollen release respectively pollen production is more sensitive to meteorological factors than the allergenicity [24][25][26]. Nevertheless for ragweed it could be shown that under elevated drought stress the expressed sequence tags (ESTs) encoding allergenic ragweed proteins increased, thus allergen content tended to increase as well [27]. Another study on Arabidopsis and rice revealed that pollen allergens tended to be part of metabolic processes in the pollen cell wall and part of stress responses [28]. These latter two studies indicate that grass pollen allergens might be enhanced under stress conditions as well.
In general, up to 95% of patients allergic to grass pollen possess IgE specific for group 1 allergens and 80% for group 5 allergens, thus these two groups make up the major grass pollen allergens [29][30][31]. Comparing the analysis of group 1 allergens (Phl p1) and group 5 allergens (Phl p5), the allergen quantification is much easier for group 5 allergens. Phl p1 reaches high homology in various grass species, but the immunodominant positions of the amino acids are different. In consequence the immune response to group 1 allergens might differ between grass species which makes an investigation of Phl p1 quite difficult [32].
The impact of drought and elevated temperature primary on the allergenic potential of different grass species and respective cultivars has up till now only been little examined. In this context it still needs to be clarified whether plant stress induces higher allergenic potentials in grass pollen. To quantify the impact of drought and elevated air temperature, this study conducted a one-year field experiment focusing on the effects of warming and drought on the phenological development, pollen weight, protein content and group 5 allergen content (Phl p5) of selected cultivars from the grass species timothy and perennial ryegrass. We hypothesize that dry conditions and elevated temperatures hamper pollen development and increase the allergen content due to plant stress.

Investigated grass species and cultivars
The following grass species and associated cultivars were selected: perennial ryegrass (Lolium perenne L.): Honroso, Borsato, Indra and Ivana; timothy (Phleum pratense L.): Comer, Lischka, Classic and the timothy grass mixture solely from the provenance Giggenhausen (48.363239˚N, 11.649388˚E); and cocksfoot (Dactylis glomerata L.): Musketier, Revolin, Diceros and Lidaglo. The timothy grasses from the provenance Giggenhausen are naturallyoccurring grasses which are summarized in the following as one group under the name of their provenance. Due to the use of relatively old seedling material in the Lidaglo cultivar with lower germination rates and slower development in the Revolin and Diceros cultivars, the pollen production of cocksfoot was insufficient, whereby appropriate pollen amounts were produced only on the control plots by the cocksfoot cultivar Musketier. Therefore cocksfoot was excluded from the analysis since no meaningful comparisons were feasible. In the beginning of the survey, the cultivar Classic (timothy) developed more slowly in its early growing stage, i.e. seemed to be undersized and less vital, therefore it was excluded from the phenological, height and photographic recordings. Later on, Classic recovered in all treatments so that this cultivar was included into the pollen sampling and the following analytics.

Experimental site
The experimental study site (Fig 1) is located at the research station Dürnast (48.404457˚N, 11.690464˚E; 445 m a.s.l.), 50 km north east of Munich, southern Germany. The site is part of the well-maintained Gewächshauslaborzentrum Dürnast of the Technical University of Munich. In total the experimental area covered 80 m 2 in which 144 plots (each 65 x 50 cm) were installed, 36 for each of the four treatments (Fig 1 and Table 1). The original loamy soil of previous experiments [33] had been exchanged by more sandy soil material with higher drainage capacity [34]. In the beginning of 2014, this material was classified as loamier sand with a 70% proportion of sand and low phosphorus and potassium content. Before sowing of the grass cultivars in autumn 2016, 5 cm of humus was applied. For each plot 0.54 g pure seeding material from one cultivar was mixed with soya grist and then equally distributed over the respective 65 x 50 cm area. Each treatment comprised respectively four different cultivars each from the three grass species timothy, perennial ryegrass and cocksfoot (4 � 3 = 12 plots). Within each treatment three repetitions were conducted on a total of 36 plots (see Table 1).

Treatments
During the investigation period from May 19 until June 23, 2017, the simulation of drought (-stress) and/or elevated air temperature was conducted in three different treatments and one control (Fig 1). Besides the control, the experiment comprised the treatments drought (rain exclusion by a rainout shelter), warming (elevated air temperature) and drought + warming (rain exclusion and elevated air temperature combined).
The elevation of air temperature was achieved by a micro-capillary warm water system (type P.VS30, Beka Heiz-und Kühlmatten GmbH, Berlin, Germany). The capillary mats (5,700 mm length, 630 mm width, capillary tube diameter 4.5 mm, distance between tubes 30 mm) were fixed on a wooden frame in 20 cm height. One single mat covered 9 plots in a row https://doi.org/10.1371/journal.pone.0248759.g001 (Fig 1). The warming system was first installed in May 2017, when first grasses reached the height of 20 cm, in order to prevent shading effects as long as possible. In the treatment plots warming and warming + drought, the air temperature at 20 cm height above ground was increased on average by 0.87˚C during the investigation period from May 19 until June 23, 2017.
For the treatments drought and warming + drought, a transparent mobile rainout shelter simulated drought by omitting any precipitation during the investigation period (Fig 1). It was controlled by a rain sensor operating the shelter as soon as the first rain drop hit the sensor. In turn, the shelter reopened again when no further drops hit the sensor.

Plant development, soil moisture and meteorological data
After sowing in autumn 2016 plants developed well and plots were evenly covered. During the initial growing phase in 2016 and also in spring 2017 all plots were irrigated regularly by a lawn sprinkler for respectively 30 min in the morning and afternoon in order to facilitate plant development. From mid of May till end of June 2017 vegetative and reproductive phenological microstages of each plot were recorded with the expanded BBCH code on a weekly basis following [35]. Observations included all microstages between macrostage 4 (booting), 5 (inflorescence emergence, heading), 6 (flowering, anthesis), 7 (development of fruit), 8 (ripening) Table 1. Experimental setup in Dürnast comprising respectively four different cultivars from the grass species timothy, cocksfoot and perennial ryegrass under three treatments and the control, three repetitions each, cocksfoot was excluded from this study since the pollen production was insufficient.

Control
Warming Grass species and 9 (senescence). The (micro-) stage was recorded for each plot on a weekly basis. Linear interpolation was used to receive the exact starting dates of flowering (BBCH 61). The average culm height for each plot was taken in parallel to the phenological recording for each plot. On the respective harvesting dates of the cultivars, inflorescence lengths of the cut culms were measured separately for each plot.
During the investigation period digital images of each plot were taken on May 30 and June 12, 2017. The green value (DN, digital number) of each RGB image ( � .jpeg) was extracted and analyzed using the package Fiji [36] which is based on ImageJ [37]. For the interpretation, green values were regarded as proxy for plant vitality.
Between May 19 and June 23, 2017, soil moisture was recorded two to three times per week in depth levels of 100, 200 and 300 mm with a soil moisture sensor (PR2 /6 SDI-12, HH2 Moisture Meter, Delta-T Devices, Cambridge, UK) at 36 spots equally distributed among the treatments where measuring tubes had been embedded in the soil. Due to very low soil water content and clear signs of dehydration of the grasses, all plots were irrigated with a watering can by 1.6 l per plot on June 01, and by 3 l per plot on June 09/14/21, 2017, respectively.
Meteorological data (air temperature, air humidity, and precipitation) in hourly resolution to characterize the growing conditions were obtained from a nearby climate station (Weihenstephan-Dürnast, location 48.4029˚N; 11.7305˚E, distance to field site 385 m) of the German Meteorological Service (DWD). In addition, air temperature and relative air humidity in 20 cm above ground were directly measured by 12 sensors at the site, equally distributed among the treatments.

Collection of grass pollen and pollen count per blade
In-situ pollen collections are very likely to be influenced by humidity (e.g. caused by rain events), mildew and insects after plants have been covered with collective containers [38]. To overcome these issues, grasses were kept in climate chambers during the actual pollen release. Accordingly, when first flowering (BBCH 61) was observed for a cultivar/treatment, bunches of 10-20 individual plants per plot were harvested and inflorescences were covered with pergamin bags and closed at the bottom. After attaining the full flowering stage (BBCH 65) under fixed conditions in a climate chamber (day/night cycle: 14 hours day at 23˚C, 10 hours night at 15˚C; air humidity 45%), pollen was extracted from the pergamin bags by shaking and subsequent removal of anthers and culms. Additionally, the number of culms per bag was counted and the total pollen weight was determined by an electronic balance (XS204DR, Mettler Toledo GmbH, Gießen, Germany). This method for pollen collection and isolation [39] was consistently applied for all samples. However, it cannot be excluded that single pollen still adhered to the bags or were not released by the anthers. To preserve pollen in the same fresh conditions before analytical testing, they were stored at -20˚C.

Extraction and determination of protein content
Following Jung et al. [38], pollen grains were mechanically extracted. Total soluble protein content was quantified using BCA test [38]. The reagents (BCA solution and copper sulfate) were ordered from Sigma (B9643; C2284) and for the standard curve Albumin from Serva (11930) was used.

Grass pollen weight and allergen quantification
Grass pollen weight was measured by dissolving 5 mg pollen grains of each grass cultivar sample immediately before the measurement in 250 μl PBS from which in turn 10μl were counted (n = 4) using an automated cell counter (TC-10, Bio-Rad Laboratories GmbH, München, Germany) [40]. All samples were counted within one day.
For the quantification of group 5 allergen content a sandwich ELISA (Allergopharma GmbH, Reinbeck/Hamburg, Germany) [41][42][43] was used, with a sensitivity of 1 ng/ml and precision of ±10% [38,44]. A standard curve was set up by the timothy (Phleum pretense) group 5 allergen Phl p5 (Allergopharma GmbH) covering a concentration range of 1 to 1000 ng/ml [38,45]. The epitopes present on the grass pollen allergens Phl p5a and Phl p5b were fixated with the monoclonal antibodies MoAb 1D11 and MoAb B01 (Allergopharma GmbH) and spectrophotometrically visualized with a chromogen present on the biotinylated MoAb B01. Since the group 5 allergens are homologous proteins in grass pollen, the same antibodies could be used for all species of Poaceae [46,47].

Statistical analyses
For the parameters air temperature at 20 cm height, soil moisture, green value, pollen production, pollen weight, protein content and allergen content the Shapiro-Wilk normality test was performed. Since the parameters were not normally distributed (p-value of Shapiro-Wilk test < 0.05), non-parametric tests were used. In order to compare more than two groups, the Kruskal-Wallis test was used to check for significant differences (p-value < 0.05) between the groups (e.g. temperature at 20 cm height in four treatments). Afterwards, the pairwise Wilcoxon test was chosen to identify significant differences between single pairs of more than two groups. Due to the limited number of samples for the respective treatment and cultivar, differences between treatments were tested on the species level, and not separately for each cultivar. Correlation analysis for non-parametric data was calculated by Spearman's Rank Correlation. P values smaller than 0.05 were considered to be statistically significant. All data were analyzed with R [48] using the packages ggplot2, latticeExtra, tidyr and ggpubr.

Growth conditions
According to the data recorded by DWD between April 24 and July 02, 2017, the average air temperature was 15˚C and the average precipitation per day was 2.8 mm (Fig 2). There was a longer dry period between mid-May until end of May with less than 1 mm precipitation in total.
On the plots with warming system air temperature in 20 cm height above ground was on average 0.87˚C higher and relative air humidity 1.34% lower compared to the treatments without warming between May 17 and July 02, 2017 (Fig 2), but these differences between the treatments were not significant (p = 0.325).
Soil moisture in 100, 200 and 300 mm depth was on average 28% lower on the plots where rain was excluded by the rainout shelter (p<0.001). During the months May and June soil moisture in 100 mm depth was on average 0.07 m 3 m -3 in the warming treatment, 0.05 m 3 m -3 in the drought treatment, 0.06 m 3 m -3 in the warming + drought treatment and 0.08 m 3 m -3 in the control treatment (Fig 2). The minimum soil moisture of 0.04 m 3 m -3 was reached on June 08 for the drought treatment. Afterwards, around June 14 there was a peak in soil moisture for all treatments due to the irrigation on June 09/14/21, 2017 (see section 2.4).
Except Ivana, all other cultivars of perennial ryegrass tended to have a slightly faster phenological development in the treatments warming, and drought + warming (p = 0.47). In general, the cultivars of timothy started flowering later than perennial ryegrass. Among timothy, Giggenhausen was the first starting to flower on June 13 (DOY 164), second the cultivar Lischka (between DOY 164 and 167), third Comer (DOY 167) and latest Classic (between DOY 167 and 171) (Fig 3). Except for the drought treatment, phenological curve progression was mostly similarly among the other treatments (p = 0.77).
Until end of May the average height did not differ between perennial ryegrass (28.5 cm) and timothy (27.6 cm) (p = 0.74) (Fig 4). In the period May 17 to June 23, perennial ryegrass grew on average in height by 17.2 cm and timothy by 16.7 cm. Ivana (perennial ryegrass cultivar) had the largest total height (47.8 cm) among all investigated cultivars, whereas Borsato had the lowest total height (33.2 cm). There were no significant differences (p = 0.75) among the treatments w.r.t. height growth between May 17 to June 23: the control treatment had on average the strongest increase in height (19.3 cm), followed by the warming (16.6 cm), drought (16.2 cm) and warming + drought (15.2 cm), respectively. Up and downward fluctuations in the height progression can be explained by the (partial) removal of plant individuals at harvesting.
The inflorescence and spikelet length of perennial ryegrass was on average 9.1 cm and for timothy 9.6 cm. There were no significant differences between the treatments. When comparing the treatments, the largest growth in height (37.8 cm) and inflorescence length (13.5 cm) was reached under the control treatment for perennial ryegrass whereas timothy had the largest height (43.5 cm) under the warming treatment and highest inflorescence length under the drought treatment (4.0 cm). There were no significant differences between the treatments for the inflorescence length.

Fig 2. Meteorological parameters during the investigation period (May 17 to July 02, 2017), top to bottom: Soil moisture [m 3 m -3 ] in 100 mm depth (measuring points were only connected for consecutive days, error bars indicate the standard error of 7-9 measurements each), mean daily temperature [˚C] as deviation from the average of all four treatments at 20 cm height (colors represent the treatments, error bars indicate the standard error of 3 measurements each), daily sum of precipitation [mm] and irrigation [mm] (in light blue), daily air temperature [˚C] obtained from the German Meteorological Service (DWD) (in red).
https://doi.org/10.1371/journal.pone.0248759.g002

Pollen analytics
Weight per grain. On average the pollen weight per grain was 30% higher in timothy (18.3 ng) compared to perennial ryegrass (14.1 ng, p = <0.01) (Fig 5). Among the observed timothy cultivars, Lischka had the highest pollen weight (20.1 ng) and Classic the lowest (16.6 ng). Among the perennial ryegrass cultivars, Ivana had the highest (19.3 ng) and Indra the   lowest weight per grain (10.8 ng). There were significant differences between the treatments for the perennial ryegrass cultivars (p<0.05) and for timothy cultivars (p<0.01).
Based on the limited number of samples per treatment and cultivar, statistical tests comparing treatments were carried out for all cultivars of one species combined and not separated by cultivar. For perennial ryegrass significant differences in the weight per pollen grain could particularly be seen between the control (17.2 ng) and the treatments drought (11.6 ng) (p<0.01), and drought + warming (12.6 ng) (p<0.05). For timothy significant differences were observed between the control (21.5 ng) and the drought treatment (15.3 ng) (p<0.01).
Protein and allergen content. The average protein amount was 32% higher in timothy (2.7 ng higher per grain) than in perennial ryegrass (2.0 ng per grain, p = <0.01) (Fig 5). Regardless of the species, Indra had the lowest protein content (1.5 ng) and Comer the highest (2.9 ng). Protein amount for the timothy cultivars significantly differed between the treatments control (3.1 ng) and drought (2.4 ng) (p<0.01). For the perennial ryegrass cultivars there was no significant effect of the treatments.
The absolute allergen content was on average more than three times higher in timothy (0.032 ng) compared to perennial ryegrass (0.010 ng) (Fig 5). The perennial ryegrass cultivar Hornroso had on average the lowest allergen content (0.008 ng) and the timothy cultivar Lischka (0.035) the highest.
The absolute allergen content (ng) was consistently and significantly higher in the control (0.040 ng) in comparison with all treatments for the timothy cultivars (average 0.029) (p<0.01). In the case of the perennial ryegrass there was a significant difference between the treatments drought (0.008 ng) and warming + drought (0.012 ng) (p<0.05).
Pollen per culm. Highest pollen production was observed for the timothy cultivar Classic (on average 3.26 mg/culm) and lowest for the perennial ryegrass cultivar Indra (on average 0.11 mg/culm) (Fig 5). Timothy grass produced on average much more pollen (1.44 mg/culm) than perennial ryegrass (0.24 mg/culm, p = <0.001).

Green values
The green value of each plot was determined twice (May 30 and June 12, 2017) by image analysis. For the drought and warming + drought plots, the green value decreased on average by 10% and 5% respectively within these two weeks (Fig 6). In contrast the green values of the control and warming plots hardly changed, in case of the timothy cultivars Comer and Lischka the green value even slightly increased (Fig 6). Green values were generally lower (-13%) on the warming treatment plots due to the capillary mats (see Fig 1). The green values for perennial ryegrass significantly decreased between May 30 and June 12 for the control (p<0.001), the drought treatment (p<0.001) and for the drought + warming treatment (p<0.05). For timothy, a significant difference between green values of May 30 and June 12 could only be seen for the drought treatment (p<0.01).

Relationship between pollen characteristics and drought stress
The correlation analysis by Spearman between the soil moisture at all depths and the height of grasses revealed a significant positive correlation (p<0.05), same for the green value on June 12 (p<0.05) and the pollen protein content (p<0.05) (Fig 7). The weight per pollen grain and the protein content were highly significantly and positively correlated (r = 0.91, p = <0.001), as well as the weight per pollen grain and allergen content (r = 0.61, p<0.001), and the weight per pollen grain and the height (r = 0.33, p<0.01) (Fig 7). There was a significant correlation between pollen production per culm and allergen percentage (r = 0.46, p = <0.001) and respectively allergen content (r = 0.44, p<0.001). No significant correlation was found between the green values of May 30 and allergen content (ng) (r = 0.21, p = 0.06) whereas a significant correlation was registered for June 12 (r = 0.42, p<0.001).

Comparison with previous studies
Our results demonstrated that climate change related drought and warming clearly influence the development of grass in different ways e.g. by altering the start of flowering, plant development and pollen characteristics. At the same time, the recorded effects of drought and warming were strongly dependent on the grass species/cultivar. First of all, the results of this study are in close accordance with comparable studies. The collected phenological data are similar to previous observations on plant development and the onset of flowering [38]. Although the grass individuals observed in this study are in the juvenile state, the determined ranges for pollen weight, protein-and allergen content correspond well to those of other studies [38,42]. Solely the error bars for the pollen production are relatively large due to the limited number of samples. Our results reveal significant differences between the grass species perennial ryegrass and timothy, particularly in protein and allergen contents, which are higher for timothy than for perennial ryegrass, regardless of the treatment. This result for the protein content is consistent with a previous study [49], i.e. the intraspecific variation in pollen production and allergen content is greater in perennial ryegrass than in timothy.

Grass species and -cultivar dependent response to drought and warming
Elevated temperatures slightly advance the start of flowering of grass species, perennial ryegrass and timothy. This finding supports previous studies, where elevated air temperatures also accelerated the onset of flowering similarly [2,50]. In contrast, the plant development is slightly delayed under the drought treatment with ambient temperatures [51]. The effect of drought becomes evident in the reduction of plant vitality and the markedly change in pollen development.
Pollen weight, protein-and allergen content of the timothy cultivars in drought and warming + drought treatments were in most cases significantly lower than the respective values of the control. Since the probes for pollen weight, protein and allergen content were sampled independently using different analytical methods, a systematic impact of drought on those parameters can be concluded. This systematic effect was found exclusively for timothy, but hardly for the cultivars of perennial ryegrass. The effect itself can be explained by the steadily decreasing plant vitality, which even required additional watering of all plots in mid-June. As most of the perennial ryegrass cultivars started flowering slightly earlier and in the case of Ivana clearly earlier than the timothy cultivars, pollen development was completed before plant vitality was seriously affected. This may explain why systematic effects have been observed for timothy but not for perennial ryegrass. The decrease in plant vitality was also seen in the change of green values between May 30 and June 12 on the plots where rain was excluded. During this period, the green values fell sharply and a slightly slower growth rate was recorded, clearly a sign of water shortage impacting plant development. The comparably high sensitivity of Poaceae to drought has already been shown by [20].
The response of the perennial ryegrass cultivars to warming/drought was in comparison to timothy more dependent on the chosen cultivar, but nevertheless a tendency towards higher allergen levels were observed under all three warming/drought treatments. For other pollen characteristics, such as pollen weight, protein and allergen content, no systematic response to the treatments was observed for the perennial ryegrass cultivars. The different reactions of the perennial ryegrass cultivars could be explained by their variable environmental constraints for growth/plant development. In consequence, their individual tolerance levels to drought or elevated temperatures could be very different. This would also explain why drought led to an increase in allergenicity in most of the perennial ryegrass cultivars, whereas a decrease was observed for most of the timothy cultivars. Specifically, the allergen proportion increased under elevated air temperatures/drought for the ryegrass cultivar Indra, Ivana and Hornroso and for the timothy grass cultivar Classic, while a decrease was observed for all other cultivars. To some extent, the influence of elevated air temperature on plant development was certainly overlaid by the interspecific variability between species and cultivars. Regarding pollen production, our results suggest that under elevated air temperatures/drought there is a high dependency on the respective cultivar.
From the climate change perspective our results propose that the allergenic risk due to timothy grass pollen will reduce in the future as the absolute allergen amount decreased under the drought treatment for all four timothy cultivars and under the warming + drought treatment for three cultivars.
The response of perennial ryegrass to drought/warming depends on the cultivar; both increases and decreases of the absolute allergen amount are possible. In comparison to timothy, perennial ryegrass has a much higher number of cultivars with 246 cultivars listed by the Bundessortenamt Germany [52,53], whereby each cultivar is highly specialized to certain environmental conditions. It has to be kept in mind that timothy had up to 5 times higher allergenic content and up to 30 times higher pollen production compared to perennial ryegrass.

Suggestions for experimental setup
Seeds for our study were obtained from provenances in Bavaria. It cannot be ruled out that specimens of the same species/cultivar but from different provenances show genetic or phenotypic variations and therefore their response to respective treatments might be different. In addition, our study is limited to a one-year field experiment and therefore can only show a portion of the perennial plant's life cycle. It could be possible that the cultivar-/ species-specific patterns shown vary within the life cycle. For example, most of the cocksfoot cultivars cultivated did not flower in the first year regardless of treatment. For validation purposes the experimental setup with the same plants was repeated in 2018. Unfortunately, due to severe technical malfunctions of the climate chambers, results from 2018 were unusable for further analysis.

Projected changes in future grasslands
Past experimental manipulations of water availability in tallgrass prairie have concluded that water limitations may be important in some years but not in others [23]. Therefore more replications of this experimental approach should be performed in several consecutive years. Due to technical prerequisites and limitations our warming treatment of 0.87˚C in 20 cm height was moderate (comparable to RCP2.6), in future stronger warming scenarios as RCP 4.5 or RCP 8.5 should be considered. The generated drought in this study simulated conditions which can be already found nowadays, also within certain areas in Germany (e.g. <30 mm precipitation in May 2017 for Brandenburg and Saxony, Source: German Meteorological Service). Dependent on the soil type and if the generated drought is stronger than in this study (<0.04 m 3 m -3 ), the permanent wilting point (pF < 4.2) will be reached in certain periods, and grasses in the actual swards might not be able to develop inflorescences and emit pollen any more. Under changing climate with decreases in precipitation, the natural grass species composition will change and adapted cultivars for agricultural production have to be selected. Within the present study we show that depending on the species/cultivars a small increase in temperature and/or drought already has a high impact on the pollen allergenicity. The differences between the species-/ cultivars which revealed in some cases even opposite behavior towards drought/warming, point out the importance of grass species/cultivars selection for people allergic to pollen. In areas where e.g. fodder crops are produced agricultural aspects such as yield, quality parameters and resistance will always be prioritized, but in the development of new cultivars and release of recommendations for adjusted seed mixtures, allergenicity should be at least considered. On the other hand, in areas with extensive use such as landscape lawn or technical grassland the allergenicity of the selected cultivars should be the key factor for decision-making.

Conclusions
In this study the impact of drought and warming on pollen characteristics such as allergen amount of the two major grass species perennial ryegrass and timothy was examined.
It could be shown that the response to drought and warming is highly dependent on the species and respective cultivars, whereby both increases and decreases of the absolute allergen amount are possible, which can be explained by the cultivar specific requirements for growth conditions and tolerance to drought and heat. For the species timothy the drought and the warming treatment led to significantly lower values for pollen weight, protein content and allergen amount. In comparison, the response to drought and warming of the species perennial ryegrass was highly cultivar-specific. Based on these results, existing knowledge about different grass species and cultivars should be expanded to include the effects of drought and warming on pollen-specific traits such as allergenicity. In the long-run, the outcome of this study can contribute to the development of climate-change adapted seed mixtures which at the same time may not increase the allergenic burden. Under changing climate those aspects have to be well studied within pollen research as they are one of the key factors in public health.